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ABSTRACT 

oo 

Context. Brown dwarf companions to stars at separations of a few AU or less are rare objects, and none have been found so far around 
early-type M dwarfs (M0V-M5V). With GJ 1046 (M2.5V), a strong candidate for such a system with a separation of 0.42 AU is 
presented. 

^ Aims. We aim at constraining the mass of the companion in order to decide whether it is a brown dwarf or a low-mass star. 

I Methods. We employed precision RV measurements to determine the orbital parameters and the minimum companion mass. We then 

derived an upper limit to the companion mass from the lack of disturbances of the RV measurements by a secondary spectrum. An 
even tighter upper limit is subsequently established by combining the RV-derived orbital parameters with the recent new version of 
C^l the Hipparcos Intermediate Astrometric Data. 

Results. For the mass of the companion, we derive m > 26.9 Mj up from the RV data. Based on the RV data alone, the probability that 
the companion exceeds the stellar mass threshold is just 6.2%. The absence of effects from the secondary spectrum lets us constrain 
the companion mass to m < 229 Mj up . The combination of RV and Hipparcos data yields a 3<r upper mass limit to the companion 
mass of 112 M Jup with a formal optimum value at m = 47.2 M Jup . From the combination of RV and astrometric data, the chance 
probability that the companion is a star is 2.9%. 

Conclusions. We have found a low-mass, close companion to an early-type M dwarf. While the most likely interpretation of this 
5^ object is that it is a brown dwarf, a low-mass stellar companion is not fully excluded. 

Key words. Stars: low-mass, brown dwarfs - Binaries: spectroscopic - Stars: individual: GJ1046 - Astrometry 
1 . Introduction where M and m are, respectively, the mass of the star and com- 
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panion, and P, e and i are the orbital period, eccentricity and 



CO The paucity of brown dwarf companions to solar-like stars at inclination 

separations of a few AU or less (a canonical value of < 5AU is 

^ usually quoted), was already noted by Campbell et al. (1988) As can be seen from equation (1) the chances to detect a 

in their early precision radial velocity survey. This "brown companion object via RVs increase with shorter period (and 

m dwarf desert" is currently not well understood (see Grether & shorter separation) and higher companion mass. For example, 

O Lineweaver 2006 for an overview). Two distinctive formation the RV semi-amplitude of the stellar reflex motion caused by a 

OO mechanisms seem to be at work for planetary (M < 13 M Jup ) and 20 M Ju P brown dwarf at 1 AU from a solar-mass star is 565 ms , 

O stellar (M > 0.08 M ) companions with relatively little overlap lf the orblt 18 Clrcular and seen edge-on (i = 90°). This is two or- 

> between the two. At wide separations no "brown dwarf desert" is de rs of magnitude greater than the current state-of-the-art RV 

observed. While the frequency of brown dwarf companions sep- measurement precision of a few ms . The detectabihty also m- 

arated from their host star by < 3 AU is about 0.5%, it is at least a creas es for lower-mass stars; the same 20 M Jup at 1 AU from an 

H factor of 10 higher for separations > 1000 AU (Gizis et al. 2001 ; 0.3 Mq star would produce an RV semi-amplitude of 1009 ms 1 . 

also Neuhauser & Guenther 2004). The fact that close-in brown And in fact there is observational evidence of the existence of 

dwarf companions are rare is highly significant since the com- brown dwarf companions to low-mass stars such as the proto- 

monly employed radial velocity (RV) method to search for sub- type brown dwarf companion GJ 229B that orbits an Ml V star 

stellar companions to stars is very sensitive to such objects. The at a wide projected separation of 44 AU (Nakajima et al. 1995). 

RV semi-amplitude K of the primary (the companion is usually Few brown dwarf companion candidates in the separation 

not visible in the spectrum) is given by regime up to a few AU are known. The first such candidate was 

HD 114762 (Latham et al. 1989), and a more recent example 

K = (2nG/P) 1/3 msinz is HD 137510 (Endl et al. 2004). The masses of these candi- 

(1 - e 2 ) l/2 (M + m) 2/3 ' date objects have often not been well-determined since the RV 

method just yields a minimum mass and the astrometric preci- 

Send offprint requests to: M. Kurster sion of the available Hipparcos data (ESA 1997) is mostly not 

* Based on observations collected at the European Southern sufficient to confirm brown dwarfs and exclude stellar compan- 

Observatory, Paranal, Chile, programmes 173.C-0606 and 078.C-0829 ions (e.g. Pourbaix 2001 ; Pourbaix & Arenou 2001). Among the 
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Fig. 1. RV time series of our UVES RV data for GJ 1046. The solid line corresponds to the Keplerian solution with a period of 
169 d and with^ 2 = 12.7, 8 degrees of freedom (DoF), and p(x 2 ) = 0.123. The rms scatter of the measurements around the orbital 
solution is 3.56 ms -1 which compares well with the average measurement error of 3.63 ms -1 which is much smaller than the plot 
symbols. For comparison the second best solution with a period of 338 d is shown as a dashed line. While being the second best, 
the latter solution is clearly excluded due to its extremely large x 2 value of 10 121. 



best established brown dwarfs are the companions to the G4IV 
star HD 38529 and to the G6IV star HD 168443 with companion 
masses of 37 M Jup and 34 M Jup , orbital periods of 2174.3 d and 
1770 d, and separations of 3.68 AU and 2.87 AU, respectively. 
These masses were determined by Reffert & Quirrenbach (2006) 
who derived new astrometric solutions from the Hipparcos mea- 
surements given the precisely known RV-derived orbital parame- 
ters, i.e. period, time of periastron passage, eccentricity, and lon- 
gitude of periastron. Another example determined in a similar 
fashion by Zucker & Mazeh (2000) is the G5IV star HD 10697 
with a companion in a 1078 d orbit at a separation of 2.12 AU 
and a mass of 40 M Jup . An example for an object with a mini- 
mum mass of 9.3 M Jup that turned out to be a star with a mass of 
142 M Jup is the companion to HD 33636 (Bean et al. (2007). 

The present paper presents a new candidate for a brown- 
dwarf companion which we have found in our precision RV 
survey carried out with the UVES spectrograph at the ESO 
VLT in search for planetary and substellar companions to M 
dwarfs (see Kiirster et al. 2003). The host star, GJ 1046 (M2.5V; 
V = 11.62 mag), has no entry in the Double and Multiple 
Systems Annex of the Hipparcos data base. Comparing its V 
band and J,H,K band colours (from the 2MASS catalogue; 
Skrutskie et al. 2006) with the mass-luminosity relationships 
by Delfosse et al. (2000) there is no indication of near-infrared 
emission in excess of the scatter found in these relations. 

If the companion to GJ 1046 turns out to be a brown dwarf, 
then the system would be unique in that it would contain the first 
close-in (< 5 AU) brown dwarf companion to a main-sequence 
star of spectral type early-M (M0V-M5 V). Since the mass ratios 
of binary systems with low-mass primaries tend towards unity, 
brown dwarf companions to late-M dwarfs (>M6V) are rela- 
tively frequent (e.g. Montagnier et al. 2006). But even in these 
systems separations < 5 AU are usually found only among bina- 



Table 1. Differential RV time series measurements of GJ 1046. 



Date 




a) BJD - 


DRV 


RV-error 






2 450 000 


[ms 1 ] 


[ms- 1 ] 


2004-10 


03 


3281.83304 


-1132.0 


3.2 


2004-11 


11 


3320.62631 


411.5 


3.5 


2005-07 


27 


3578.91097 


-1670.6 


4.8 


2005-08 


26 


3608.77350 


-1524.6 


3.7 


2005-09 


11 


3624.71807 


-928.7 


4.3 


2005-09 


19 


3632.72200 


-615.2 


3.5 


2005-10 


17 


3660.68277 


492.3 


3.5 


2006-01 


15 


3750.60256 


-1787.7 


3.5 


2006-09 


15 


3993.84856 


326.4 


3.3 


2006-09 


28 


4006.77107 


823.4 


3.9 


2006-10 


05 


4013.80030 


1085.8 


3.6 


2006-10 


06 


4014.61077 


1108.1 


3.7 


2006-11 


02 


4041.65276 


1737.5 


3.1 


2006-11 


08 


4047.65399 


1673.7 


3.3 



Note: a) Barycentrically corrected Julian Date 



ries with low-mass stellar secondaries. Counterexamples are the 
M8 star LHS 2397a with an L7.5 companion at a separation of 
2.9 AU (Freed et al. 2003) and the young M6 object at the star- 
to-brown dwarf border, Cha Ha 8, orbited by a brown dwarf at a 
separation of 1 AU (Joergens & Miiller 2007). 

2. Observations 

GJ 1046 was observed with the VLT-UT2+UVES as one of the 
targets of our precision RV survey of M dwarfs in search for 
extrasolar planets (see Kiirster et al. 2003, 2006). To attain high- 
precision RV measurements UVES was self-calibrated with its 
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Table 2. System parameters of GJ 1046 



RV-derived parameters 


Orbital period P 


168.848 


±0.030 


[d] 


Time of periastron T p 








BJD - 2 450 000 


3 225.78 


+0.32 




RV semi-amplitude K 


1830.7 


±2.2 


r 1 n 

[ms l ] 


Orbital eccentricity e 


0.2792 


±0.0015 




Longitude of periastron a> 


92.70 


±0.50 


[°] 


Mass function f(m) 


9.504 


±0.024 


[10- 5 M o ] 


^ 2 (DoF=8) 


12.7 






p(x 2 ) 


0.123 






Scatter rms 


3.56 




[ms" 1 ] 


Mean error ARV 


3.63 




[ms- 1 ] 


Inferred parameters 


Stellar mass M 


0.398 


±0.007 


[M Q ] 


Minimum companion 








mass m m i„ 


26.85 


±0.30 


[M Jup ] 


Semi-major axis of 






companion orbit a 


0.421 


±0.010 


[AU] 


Critical inclination 








(for m = 0.08M Q ) i' crit 


20.4 




n 


a) Probability of i crit , p, cril 


6.3% 


(20.4° > i > 0°) 





Parameters derived from absence of companion spectrum 



Maximum companion 

mass mX 
Minimum inclination i s ? 

mm 

a) Probability of ij^, p i sp 



229 
8.7 
1.2% 



[M Jup ] 
[°] 



(8.7° > i > 0°) 



Astrometry-derived parameters 



Ascending node £1 


97.7 


formal optimum 


[°] 


Inclination i 


125.9 


formal optimum 


[°] 


Companion mass m 


47.2 


formal optimum 


[M Jup ] 


X 2 (DoF= 202) 


329.0 






Minimum inclination 


15.6 


3<r limit 


[°] 


a) Probability of fj", p t u 


3.7% 


(15.6° > i > 0°) 




Maximum companion 








mass m^ x 


112 


3cr limit 


[M Jup ] 


b) Probability of a stellar 








companion p. 


2.9% 







Notes: a) A priori probability based on the RV derived minimum 
mass and assuming random orientation of the orbit, 
b) Probability derived from the astrometric model. 



iodine gas absorption cell operated at a temperature of 70°C. 
Image slicer #3 and an 0.3" slit were chosen yielding a resolv- 
ing power of R = 100 000 - 120 000. The central wavelength of 
600 nm was selected such that the useful spectral range contain- 
ing iodine (I2) absorption lines (500 - 600 nm) falls entirely on 
the better quality CCD of the mosaic of two 4 K x 2 K CCDs. 

For a detailed description of our data modelling approach 
employed for the determination of high precision differential ra- 
dial velocities (DRV) we refer the reader to Endl et al. (2000). 
A concise summary can also be found in Sect. 4 of Kiirster 
et al. (2003). 

A total of 14 spectra of GJ 1046 observed through the iodine 
cell were obtained in 14 nights between 03 October 2004 and 
1 1 November 2006. See Table 1 for the journal of observations. 
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Fig. 2. Histogram of the mean internal RV measurement errors 
for our sample of stars. Applying a k-ct clipping procedure with 
an iterative rejection of values exceeding the mean plus 2.6 cr 
(for a Gaussian distribution the chance probability of exceed- 
ing this value in one iteration is 0.5%) we reject all stars with 
mean RV errors in excess of 6.5 ms -1 . Of the total sample of 41 
stars only 38 are displayed, because three stars with large errors 
(79, 81, and 92 ms -1 , respectively) lie far outside the displayed 
error range. They are newly discovered double-lined spectro- 
scopic binaries with such strong contributions from the secon- 
daries that the employed single-lined spectrum model fails. Six 
stars have abnormally high mean errors between 7 and 16 ms -1 
which, in one case, is the result of very low signal-to-noise ratio, 
but in the other cases could also be due to spectral contamina- 
tion from a companion. The remaining bulk of the distribution 
(mean error < 6.5 ms 1 ) has a mean of 3.91 ms 1 and a width of 
<x = 1.02 ms- 1 . 



Individual exposure time was 900 s yielding an average S/N per 
pixel between 39 and 58 for the various spectra (the median and 
mean being 51.0 and 49.1, respectively). On average our error 
of the individual RV measurements is 3.63 ms -1 for this V = 
1 1 .62 mag object. All RV data were corrected to the solar system 
barycenter using the JPL ephemeris DE200 (Standish 1990) for 
the flux-weighted temporal midpoint of the exposure as given by 
the UVES exposuremeter. For each epoch of observation proper 
motion corrected stellar coordinates were used. On 19 November 
2004 we also obtained a triplet of exposures without the iodine 
cell (exposure time 3 x 705 s) required as a template spectrum in 
the data modelling process (cf. Endl et al. 2000). 



3. Results 

Our RV time series for GJ 1046 is listed in Table 1 and also 
shown in Fig. 1 along with the best-fit Keplerian orbit yielding 
an orbital period of P = 169 d, an eccentricity of e = 0.28, 
an RV semi-amplitude of K — 1831 ms -1 , and a mass function 
f(m) = (msim) 3 /(M + m) 2 = 9.5 x lO" 5 M G (Table 2). 

Given the unknown inclination i only a minimum to the 
companion mass can be determined corresponding to the case 
i = 90°. For this we need an estimate of the stellar mass of 
GJ 1046 which we obtain from the K-band mass-luminosity 
relationship by Delfosse et al. (2000). Taking the apparent K- 
magnitude from the 2MASS catalogue (K = 7.03 mag) and com- 
bining it with the Hipparcos parallax (71.11 mas) we find an ab- 
solute K-magnitude of 6.29 mag. The K-band mass-luminosity 
relationship then yields a stellar mass of 0.398 ± 0.007 M . 
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Table 3. Simulations for the determination of an upper limit to the companion mass. 



Mass V-band Model Mean error Excess error Chance probability from 



[M G ] 


[M Jup ] 


flux ratio 


spectrum 


[ms" 1 ] 


[ms- 1 ] 


sample comparison 


excess error 














a) 


b) 


criterion 


0.00 








none 


3.63 


0.00 


0.56 


0.61 


1 


0.16 


167 


18.30 


GJ699 


3.79 


1.07 


0.50 


0.55 


0.70 


0.17 


177 


14.85 


GJ699 


3.92 


1.49 


0.47 


0.50 


0.59 


0.18 


188 


12.26 


GJ699 


4.04 


1.77 


0.47 


0.45 


0.51 


0.19 


198 


10.27 


GJ699 


4.22 


2.16 


0.34 


0.38 


0.40 


0.20 


209 


8.69 


GJ699 


4.47 


2.61 


0.25 


0.29 


0.28 


0.20 


209 


8.69 


GJ682 


4.41 


2.51 


0.31 


0.31 


0.31 


0.21 


219 


7.43 


GJ682 


4.79 


3.12 


0.17 


0.20 


0.14 


0.22 


230 


6.41 


GJ682 


5.16 


3.67 


0.13 


0.11 





0.23 


240 


5.57 


GJ682 


5.61 


4.28 


0.094 


0.048 





0.24 


250 


4.87 


GJ682 


6.10 


4.90 


0.031 


0.016 





0.25 


261 


4.29 


GJ682 


6.70 


5.63 





0.0032 





0.26 


271 


3.79 


GJ682 


7.40 


6.45 





0.00032 






Notes: a) Fraction of stars in the distribution with larger mean errors 

b) Areal fraction of Gaussian fitted to the distribution exceeding the mean error 



We then infer a minimum companion mass of m^n = 
26.9 M Jup and, from equation (1), a semi-major axis of the com- 
panion orbit of a = 0.42 AU. In order for the true companion 
mass to exceed the stellar threshold of 0.08 M Q the orbital in- 
clination ;' would have to be < 20.4°. For a chance orientation 
of the orbit the probability that i is smaller than some angle is 
given by 

p(0 > i > 0°) = 1 - cos 6 ; (2) 

hence the chance probability to have an inclination < 20.3° is 
just 6.3% making it not very likely that the companion is a star 
(see also Table 2). 



4. Spectroscopic companion mass upper limit 

An upper limit to the mass of the companion can be determined 
from the spectroscopic data by exploiting the notion that with 
increasing mass the companion would at some point become so 
bright that it would noticeably affect the RV measurements. In 
our data modelling approach a suitable indicator for the pres- 
ence of an additional perturbing signal is the magnitude of the 
RV measurement error, because all RV data are obtained assum- 
ing that the modelled spectrum is that of a single-lined binary 
in which the light from the companion can be neglected. If this 
is not the case, then the light contribution from the companion 
manifests itself as unusually large errors of the determined RV 
values. For details of the data modelling approach and the esti- 
mation of the RV errors we refer the reader to Endl et al. (2000). 

Briefly, the spectra are subdivided into a 500 chunks of 
=s 2 A width, each of which is modelled in order to yield an inde- 
pendent RV measurement. The mean of these values and its er- 
ror are taken as the final measurements and their (internal) error 
estimates. If a second contaminating spectrum is present, the dif- 
ferent spectral chunks are affected in a non-homogeneous fash- 
ion, depending on the detailed spectral line patterns within the 
chunk, with the effect that the derived chunk RV values exhibit a 
stronger scatter and hence combine to an increased value of the 
internal RV measurement error. By way of simulations adding 
faint companion spectra to the spectra of GJ 1046 we investigate 
the following two possibilities of obtaining information on the 
companion mass upper limit from the internal RV error. 



1. A comparison of the mean internal RV measurement error 
of GJ 1046 for different added companion spectra with the 
distribution of mean observed errors for the sample of mon- 
itored stars. 

2. A comparison of the mean internal RV measurement error 
of GJ 1046 for different added companion spectra with its 
original value. 

In the case of considerable contamination the average RV 
error for the spectra of the star in question stands out from the 
distribution of RV errors for the sample of monitored stars which 
is shown in Fig. 2. If only the bulk of the distribution is consid- 
ered, i.e. stars with a mean RV error < 7 ms -1 , this distribution 
is nearly Gaussian with a mean value of 3.91 ms -1 and a width 
of <x = 1.02 ms -1 . 

The average RV measurement error for GJ 1046 is 3.63 ms -1 
with the errors of the 14 individual GJ 1046 spectra ranging from 
3.07 to 4.84 ms -1 . These values are typical of the range of signal- 
to-noise values of our spectra (between 39 and 58 for GJ 1046) 
and well inside the distribution of errors in our sample of stars 
(Fig. 3). We also note that the RV errors increase with decreasing 
signal-to-noise ratio of the GJ 1046 spectra which would not be 
the case, if the errors were dominated by a perturbing signal. 

In the performed simulations we added spectra of faint low- 
mass stars to the original GJ 1046 spectra. As summarised in 
Table 3 these simulations explored the companion mass regime 
0.16 - 0.26 M Q (first two table columns) which, according to 
the V-band mass-luminosity relation by Delfosse et al. (2000), 
corresponds to companions that are factors of 18 - 4 fainter than 
the stellar primary (third column). We have used the high-signal- 
to-noise template spectra (observed without the iodine cell) of 
two other M dwarfs of this survey, GJ 699 (Barnard's star; M4V; 
mean S/N per pixel: 300) and GJ 682 (M3.5V; mean S/N per 
pixel: 230) for masses < 0.20 Mq and > 0.20 M Q , respectively 
(fourth column). For each simulated spectrum the companion 
spectrum was shifted to the appropriate companion RV for the 
probed mass value as well as scaled in flux corresponding to the 
brightness predicted by the V-band mass-luminosity relation. 

We find that a companion with a mass > 0.254 M or 
265 M Jup and a factor of 4.1 fainter in the V-band than its pri- 
mary would increase the internal error (fifth column in Table 3) 
to > 7 ms -1 which would make this star stick out from the bulk 
of the sample (Fig. 2) indicating a contamination of the spec- 
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trum. The 7th and 8th columns of Table 3 list the chance prob- 
abilities of the obtained mean errors (fifth column) from a com- 
parison with the total sample of observed stars. 

For the second possibility of determining the mass upper 
limit of the companion we assume (conservatively) that the mean 
original RV error is entirely caused by contributions from the 
companion spectrum and not attributable to photon noise or to 
effects of instrumental nature or intrinsic to the star. We then 
search for the companion spectrum (as a function of companion 
mass and brightness) whose addition to the observed spectra in- 
troduces an additional RV error (6th column in Table 3) of the 
same magnitude, i.e. it doubles the square of the errors. With an 
original value of 3.63 ms~' we search in the simulated data for 
the companion mass and brightness that leads to a mean intrinsic 
error a factor of V2 larger, i.e. 5.13 ms _1 . (We note in passing 
that this value is in the 88.4% percentile of the distribution of the 
stellar sample truncated at 7 ms" 1 ; see Fig. 2.) The 9th column 
of Table 3 lists the chance probability of obtaining the excess 
error value listed in the 6th column. 

For this increased error value we find a companion mass of 
> 0.219 M or 229 Mj up and a primary-to-secondary V-band flux 
ratio of 6.5 (cf. Table 3). 

As the mass value derived with the criterion to double the 
square of the mean internal error is lower than the one derived 
from the comparison with the star sample, we will adopt the 
value of 229 Mj up as the spectroscopic upper limit to the mass of 
the companion to GJ 1046. This value corresponds to an orbital 
inclination of 8.7°. The probability for an inclination as small 
as (or smaller than) this value is 1.2%, again assuming random 
orientation of the orbit (see also Table 2). 



5. Companion mass upper limit from a combination 
of the RV data with Hipparcos measurements 

Even if the astrometric signature of the companion is not seen 
in the Hipparcos data, Hipparcos astrometry can yield stringent 
upper mass limits on companions detected via the radial velocity 
method. 

Using the Hipparcos parallax (71.11 mas) together with the 
orbital parameters derived from the RV measurements we can 
predict the minimum astrometric signal of the stellar reflex mo- 
tion to be 3.7 mas peak-to-peak. This corresponds to the full 
minor axis of the orbit. Since the Keplerian fit to the RV data 
only permits the determination of the projected orbit of the stel- 
lar reflex motion, the true astrometric effect could be consider- 
ably higher. For the limiting inclination of 20.4° the full minor 
axis of the stellar orbit would extend 10.6 mas on the sky. 

We have analysed the Hipparcos Intermediate Astrometric 
Data for GJ 1046 (HIP 10812) using the new reduction of the 
raw data (van Leeuwen 2007 a,b). We followed the approach de- 
scribed in Reffert & Quirrenbach (2006) by keeping those of 
the orbital parameters that are known from the analysis of the 
RVs fixed and varying only the inclination and the ascending 
node while fitting an astrometric orbit to the abscissa residu- 
als. Additional free parameters in the fit were a correction to 
the mean position, mean proper motion and parallax of the star. 
The result is shown in Fig. 3 (left panel). 

The formally best fit to the Hipparcos data is achieved with 
an inclination i = 125.9° (i— 90° = 35.9°) corresponding to a true 
companion mass of 47.2 Mj up pointing at a brown dwarf com- 



panion (Table 2). [^J However, an F-test measuring the variance 
improvement yields a probability of 17% for the detection of the 
astrometric orbit implying that is has not been detected with sig- 
nificance. This can also be seen in Fig. 3 (left panel), where the 
ascending node is completely undetermined since the 2 and 3cr 
confidence contour levels span the entire parameter range. 

In the right panel of Fig. 3, the x 2 value is shown as a func- 
tion of inclination only, together with the lcr and 3cr confidence 
regions for the inclination. The 3cr (99.73% confidence) lower 
limit to the inclination is i - 15.6° implying a 3cr upper mass 
limit for the companion of 1 12 Mj up . 

Therefore, a stellar companion cannot be fully excluded, 
even though it is unlikely. From the astrometric solution the 
chance probability for the companion to have a stellar mass, or 
equivalently, for its inclination to be either i < 20.3° or > 159.7° 
is 2.2% and 0.7%, respectively, corresponding to a combined 
chance probability of 2.9% (see also Table 2).F] 

6. Conclusions 

We have presented the discovery of a probable brown dwarf 
companion to an M dwarf with an orbital period of just under 
1/2 year and a star-companion separation of 0.42 AU. Our RV 
measurements provide a lower limit to the true companion mass 
of 26.9 Mj up and a chance probability of just 6.2% that the com- 
panion is actually a star. From the absence of any indications of 
a secondary spectrum in our data we can place an upper limit to 
the companion mass of m = 229 Mj up . 

Combining our RV measurements with the Hipparcos 
Intermediate Astrometric Data from the recent new reduction 
by van Leeuwen (2007a,b) we find a formal best-fit compan- 
ion mass value of 47.2 Mj up , but pertinent to a model that is not 
significant. However, the same data allows us to place a much 
tighter companion mass upper limit of 1 12 Mj up at 99.73% con- 
fidence. This mass upper limit still allows a stellar companion, 
but with a low probability. From the astrometric analysis the 
chance probability that the companion mass exceeds the stellar 
mass threshold is 2.9%. 

If the brown dwarf nature of this object can be fully estab- 
lished, e.g. from future astrometric measurements, it would be 
the first genuine brown dwarf desert object orbiting an early-M 
dwarf. 
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